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Abstract 

The phase diagram of an SU(2)l X SU(2)r lattice Higgs- Yukawa model with 

finite A is constructed using mean field theory. The phase diagram bears a 

superficial resemblance to that for A = oo, however as A is decreased the 

paramagnetic region shrinks in size. For small A the phase transitions remain 

second order, and no new first order transitions are seen. 
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Recent experimental evidence of a top quark with mass ~ 175 GeV [1] indicates that its 
Yukawa coupling is of order 1, possibly making nonperturbative effects significant. Lattice 
Higgs- Yukawa theories provide a way of studying non-perturbative physics of theories con- 
taining interacting scalars and fermions, although the technical problems associated with 
chiral fermions restricts us to vector theories. 

Most of the work on Higgs- Yukawa theories has dealt with the limit in which the Higgs 
field is radially frozen (A = oo). Such models have been studied both analytically f^-||| 
and using Monte Carlo || for various gauge groups. Radially active (finite A) Higgs models 
without fermions have also been thoroughly examined H0], however relatively little work 
has been done on the problem including fermions ||. In this paper I estimate the phase 
diagram of the SU (2) L x SU (2) R Higgs- Yukawa theory for the full range of A and the Yukawa 
coupling. The calculations are performed using the mean field approximation (MFA) ||, 
with the fermions included in the manner used by Stephanov and Tsypin for the radially 
frozen theory ||. 

The action for the model in d dimensions is S = Sp + Sh + Vh, 

x 

= Y,MKxy + y{P R px^ x + PlPx^DK^V 

xy 

= ^24> x M xy tp y (1) 
xy 

Sh = -kY, \HPx&xP*+^x+» + Pv+ftUnP*®*) ( 2 ) 
x,n * 

Vh = Y,v(p*) = Y,KpI-i) 2 + Px- (3) 

X X 

The fermion field ip x is an SU(2) doublet, P RL = |(1 ± 7 rf+1 ), and k is the scalar hopping 
parameter. The scalar field has been separated into its magnitude and angular degree of 
freedom, p x and Q x respectively, where $ x is a 2 x 2 SU(2) matrix. We will also use the 
0(4) field 4> x satisfying (f) x (f) x = 1. The two notations are related by $ x = (p x T k where 
T k = (1, ia). Integrating out the fermions gives the effective action for the scalar field 
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S ef f = Sh + Vh — N F In det M, (4) 

where N F has been introduced as the number of fermion species, each having the same 
action Sp- The partition function for the effective scalar theory is 

Z = I I[pldp x f HDfcexpi-Seff) (5) 

J x J X 

where D(f) x is the 0(4) invariant group measure for A '" 



x m 

rk 



The variational form of the MFA is employed by introducing a parameter H , then 
adding and subtracting the trial action p x <fi x H k to obtain 

exp(-F) = / \{pldp x f Y[D ( f> x ex V (-S eff + J2p^xH k -J2p^xH k ). (6) 

X X X X 

This gives the variational limit on the free energy F, 

F < F var = (S H - N F In det M) H + (£ H k cj) k x p x ) H - In Z H , (7) 

X 

(A)h = Z H l f Y[d Px ple W (-V(p x )) f YlD^expC^H^xPM, (8) 

X X X 

Zh= I \{dp x pleM-V{Px)) [ n^exp(£# fc ^). (9) 

J x J x x 

F var is then minimized with respect to H k . The integration over p x is necessary since for 
finite A the magnitude of the scalar is no longer constrained to be 1. At small A the shallow 
scalar potential causes fluctuations about p x = 1 which alter the results considerably. 
For small Yukawa coupling the fermionic determinant is calculated by expanding in y 

(\ndetM) H = \ndetK xy - £ V E C-'CK^-i^L^ 

n=2,4,... 71 xi...x n 

winded- £ U^f) n E l\r:,y--l\ r ,,.JM*,X-.y^ J , .'K J u ■ (10) 



n=2,4,... n\ d / Xl ... Xn 

The second line uses the approximation 

K£ = ^K ay + 0{l/*) i (11) 

which is justified since the MFA is itself only accurate to order 1/d. The determinant is then 
represented by a sum of closed hopping diagrams connecting sites associated with alternating 



$'s and $f's. Previous work has approximated the fermionic determinant by including an 
infinite subset of hopping diagrams [||] , or expanding Eq. fllCf ) to some finite order in y j| . 
We will adopt the latter approach. 

The A = oo case provides a hint in choosing how far to carry out the expansion of the 
fermionic determinant. Since the radially frozen theory in the MFA agrees quite well with 
Monte Carlo results when the free energy is expanded to order y 4 , we will evaluate Eq. 
to order y 4 . It is a tedious though straightforward exercise to enumerate all such hopping 
diagrams in d dimensions to obtain 

2<i/2 

(lndetM)^ « In det K + —y 2 tr(p x <S> x } H (p x $ x ) H - 



2<i/2 



+2(2d - l)tT(p x %p x $ x ) H {p x <$> x ) H {p x ® x )H 

+2(d-l)tr(p x $ x ) 4 H 



(12) 



where the traces are over SU (2) indices, and the indices x and y simply indicate which $'s 
are distinct group elements. 

The various quantities of the form (A)h in Eq. include group integrals, and integrals 
over p. The group integrals are calculated to second order in H because that is all that is 
required to find second order critical lines. We first compute the 0(4) group integrals by 
taking derivatives with respect to H of 



J D^eMH^Px) = 2n 2 [I (H Px ) - h(Hp x )] 



2vr 2 (l + lplH 2 + ^p x H 4 ) + Om (13) 



where H = V H k H k and I n (x) is the nth order modified Bessel function. We find 



J D<f> x eMP*H n rM = YP* Ht + °(# 3 )' ( 14 ) 
J D(f) x exp( Pl .iPC)0^ = ^n 2 plWW + fi*(± + + 0{H 4 ). (15) 



The corresponding SU(2) group integrals are 



tr($ x $ x )„ = -27r 2 + 0(# 4 ), 



71" 



(^ x ) H = t—Hp x + 0(H 3 ) 



(16) 
(17) 



In Eq. [T7| the gauge has been fixed to H l = H5 l °. 

In addition, the slightly more complicated group integral involving four $'s is needed. 
This is calculated as 

tr / D<p x D<P y exp(p x H n r x +P y H n r y )®M$M 



1 

2 



y)*ryry 



The second line relies on the fact that the trace must be separately symmetric in both i, k 
and in j, I. The sum over 0(4) indices is computed using Eq. [15] and by again choosing the 
H l = H5 M gauge. 

For computing the integrals over p it is useful to define 



Pr, 



dpp n eM-V{p)) 



-(2A) exp ----- 



(19) 

where D »+i is the parabolic cylinder function. The single site partition function is then 
given by 

POO 

Z x = 2tt 2 / dp x pleM-V(p x ))[I (p x H) - h{p x H)\ 
Jo 

= 2tt 2 / d PxP leM-V( Px ))[l + -H 2 pl\ + 0(F 4 ) 
Jo 8 



2tt 2 (P 3 + -H 2 P 5 ) + 0(H*). 

o 



(20) 



Similarly, the integrals making up Eq. [7] are given by 



tr(p x §p x $) H 



P 



2P 

Pr, P? 
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H + 0(H 3 



tr(p x <S> x p y &p x $ x p y &) H 



Ps\ 2 , jp P&Pr-SPi +Q{H 4 ) 



P, 



12Po 3 



(21) 
(22) 
(23) 



Substituting Eq.'s |2T] - p3| into Eq.'s [L2| and [7| gives F var . 

There are four phases: 1) the ferromagnetic (FM) phase with 7^ 0, 2) the param- 
agnetic (PM) phase with (<j) x ) = 0, 3) the antiferromagnetic (AFM) phase with (£, x 4> x ) 7^ 
(where £ x = ^—\Y 1+X2+ - +Xd ), and 4) the ferrimagnetic (FI) phase in which (0*) 7^ and 
(Jix'fix) 7^ 0- I n the MFA these phases are indicated by the value of H k which minimizes 
the variational free energy. In the FM phase H k 7^ 0, while for the PM phase H k = 0. To 
see the AFM phase the action must be transformed so as to make the staggered magnetiza- 
tion (£e0z) accessible. Since the action is invariant under <p x — > £ x (f> x , ip x — > exp(iii2^ x /2)ip x , 
y — > iy , k — ► — k, the action obtained by this transformation will have AFM order if H k 7^ 0. 
Therefore, the MFA variational action for the AFM is given by Eq. [12] with y — > iy, k — ^ —k. 
The FI phase is indicated by the existence of FM and AFM order at the same point. 

The strong coupling regime is reached by expanding the fermionic determinant in 1/y 
rather than y. In this expansion K xy appears rather than K~y so the 1/d approximation of 
the propagator is not needed. Aside from this difference, the free energy is computed in the 
same manner as the weak coupling version. Therefore the strong coupling variational free 
energy is obtained making the replacement y — > d/2y in Eq. |T2|. 

The second order phase transitions are found by setting ^pF var \H=o = 0, which is why 
it was sufficient to compute F var to 0(H 2 ). For weak coupling this gives 

FM - PM : Ko = A _ ^ILJ + (L - iV , (24) 

P 5 d d 2 y P 3 d 3 \2d ) y y ' 

AFM - PM : Kc = - A - 2 ^Ly> - ( i- - lV. (25) 

P b d d 2 y P 3 d 3 \2d ) y v ' 

The figures illustrate the above results. Fig. 1 shows function of A at y = for 

the FM-PM transition. (At y = 0, k c for the AFM-PM transition is simply the negative 

of k c for the FM-PM transition.) In agreement with previous results PJlQfl, n c —>■ 1/8 as 

A —>■ 0, and k c —>■ 1/4 as A —>■ oo. Fig. 2 shows the complete phase diagram for A = 1, which 

appears qualitatively similar to the A = oo case. The one difference from the A = oo case is 

that the PM region has shrunk. Fig. 3 shows the Kj C S cLS cL function of y for several values 

of A between 0.01 and 10. As A is decreased, the phase diagram remains qualitatively the 
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same as the A = oo case, although the PM region shrinks in both the y and k directions. 
All phase transitions remain second order, and there is no evidence of new phases. Also, k c 
saturates for A = 0.01, 10, with little change in k c as A — > 0, oo. 

In addition to the MFA errors of order 1/d 2 , the errors in the expansion of the fermionic 
determinant are of order y 6 . For y ~ 1 the results should break down, although comparison 
of the A = oo results with those from Monte Carlo show good agreement through the 
intermediate coupling region. 
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FIGURES 



FIG. 1. k c for the FM-PM transition at y = 0. 

FIG. 2. The phase diagram for A = 1. All phase transitions are second order. 

FIG. 3. Second order k c 's for various values of A as indicated by the labels. The A = 0.01, 10 
lines are very nearly the same as those for A = 0, oo respectively. 



9 



